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The performance of AlGaN-GaN high-electron mobility transistors 
(HEMTs) is influenced by the self-heating phenomenon, which leads to the 
power dissipation that is related to the increase of the local temperature of 
the device. The study of this increasing of the temperature is executed under 
different parameters, namely, low field mobility, velocity saturation, thermal 


conductivity of the substrate. A thermal model is developed to study the 


effect of this phenomenon on the current-voltage characteristics. Among the 
Keywords: techniques to minimize the increase in the local temperature, we based on 
the good choice of the substrate used in the transistor. To highlight this 
2-DEG j i 
proposal model, we have made a comparable study between the substrates of 
AlmGal-mN/ GaN HEMTs silicon and sapphire. Our analytical results are in a good agreement with 
Analytical thermal model 


published experimental data. 
Power dissipation 


Self-heating This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Azzeddine Farti 

Laboratory of Materials Physics, Microelectronics, Automatic and Thermic 
Faculty of Sciences Ain Chock Km 8, University Hassan II 

Casablanca 5366, Morocco 

Email: azzeddine.farti-etu @etu.univh2c.ma 


1. INTRODUCTION 

The semiconductor devices has contributed in the development of information technologies and 
communication. An excellent property of electronic devices based on semiconductor materials of the group 
II-V, among them a high electronic mobility and direct bandgap. In addition, it’s able to operate at high 
frequencies [1]. The high-electron mobility transistors (HEMT) based on AlGaN-GaN have become very 
strong in recent years for the operating at high temperature, high frequency and high voltage [2]-[6]. This is 
due to the good properties of gallium nitrides (GaN), among them; we note that it has a wide band gap, a 
higher mobility of electron, as well as a very high velocity of saturation [7]. 

However, the operation of high voltage in high performance GaN integrated in circuit for example, 
low-noise amplifiers voltage controlled oscillators have been indicated [8]. This is due to a phenomenon that 
we can call it the self-heating. The electrical properties of the transistor HEMT can affect too such as, 
velocity of saturation, electron mobility, pinch voltage, breakdown voltage and band gap [9], [10], that is not 
only influence the electrical characteristics negatively, but also degrade the lifetime of the device according 
to the Arrhenius equation [11]. One of the factors that have a very interesting influence on the electrical 
characteristics of the transistor HEMT we mention the self-heating. That is motivates us to develop a 
theoretical thermal model, in order to study this phenomenon and to look for the most adequate substrate for 
a minimum power dissipation. In this work, we develop a thermal model, and we verify the evolution of the 
current voltage (J - V) characteristic under the influence of temperature, then we compare it with the 
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experimental results. Finally we use two substrates sapphire and silicon, and we study the variation of the 
electrical current-voltage (J- V) with a fixed the value of gate voltage applied V}s. 


2. Model 
2.1. Charge control model 

The Figure 1 presents the fundamental structure of transistors HEMTs AlGaN-GaN, that is 
considered in this analytical model. The three main materials that make up the structure of HEMT are: 
material with a low bandgap, material with a wide bandgap, and substrate. The formation of concentrations 
of two dimensional electron gas (2-DEG) carriers, at the AlGaN-GaN hetero-junction, it can be given by the 
following form [12]. 


Source 


UID —Al,,Ga;_,,N (Spacer layer) 


UID - GaN EES 


Figure 1. The fundamental structure of HEMT AlGaN-GaN 


2 
-K2 (+ [K3 (r)+4.K4(1).(Vgs—Ven(m)-K (7) 
2.K4(T) 


ns(T,m) = (1) 


With m present the molar fraction of Aluminum in AlGaN-GaN, and K,, K, and K3 are three 
different parameters depend on the temperature T. 


q.d 
e(m) 


With: K4(T) = K3(T) + (2) 
Where, Z e(m) and q are the gate width, the dielectric constant of Al,,Gaji,N and the electron charge 
respectively, (d = da + dep + dsp) da is the doped layer thickness of the n-AlGaN, d,, and dep are the 
thickness of the spacer layer and the cover layer of the undoped layer thickness of d’ AlmGaz-mN, respectively, 
Vin (m) is threshold voltage is given by [13]. 


Vin (m) = (m) — AEc(m) — #44 x (1 422) ez q 3) 


2.e(m) da e(m) 
With Na (m) AEc, are represent the the doping density of the AlGaN layer, the schottky barrier 


height, and the conduction band discontinuity between AlGaN and GaN respectively, opz(m) is the charge 
density induced by the total polarization can be calculated as (4) [14]. 


Opz (m) = [Ppp (m) + Pa eine (m) S P eas nO) (4) 


Where Pig iinĂ (m) denotes the spontaneous polarizations of AlmGaz-mN, and, Piac y (0) denotes the 
spontaneous polarizations of GaN, and Ppp (m) denotes the piezoelectric polarization defined as (5) [14]. 
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Pié ca, mum) Pour 0<m<0.38 
Ppp(m) = 4 pP? ca, „n(233— 3.5m) Pour 0.38 <m < 0.67 (5) 
0 Pour 067<m<1 


The spontaneous and piezoelectric expressions are presented in Table 1. 


Table 1. Parameters of the Al,,Gay..N-GaN heterostructure 


Description Parameter Expression 
Spontaneous polarization of AlGaN Pa ede ow (m) C/m?  —0.09m — 0.034(1 — m) + 0.019m(1 — m) 
Piezo-electric polarization of AlGaN Pie ea ot (m) C/m? mPia [Ep(m)] + (1 — m) PR [e, (m)] 
Piezo-electric polarization of AIN PEN Chm? —1.808¢e,(m) + 5.62422 (m) 
Piezo-electric polarization of GaN Pew C/m? —0.918£, (m) + 9.54123 (m) 
The basal strain €p(m) (ee a(m))) (m) 
The lattice constant a(m) A (0.31986 — 0.00891m) 


2.2. Current—voltage characteristics (Ids - Vds) 

By giving a fixe value of gate voltage V,, superior of the threshold voltage V,,, the movement of 
electrons at the interface of AlGaN-GaN is accompanied by the application of the drain voltage Vg, greater 
than zero. The generated of the drain current Ias is given by (6) [15]. 


av Kp dng(T, 
las(T,m, x) = Z.q.4(x) (ns (T, m) SE + TEn) (6) 


Where, Z, V-(x) and Kp are the gate width, the channel potential at position x and Boltzmann’s constant 
respectively, u(x) is the field dependent of the mobility determined by (7). 


u(x) = E 7 aV cx) (7) 
(aliea a dx ) 


Where, Vsat Ho Ec are indicate the saturation drift velocity, the low-field mobility, the saturation electric field 
respectively. By variation the bias voltage Vgs, the HEMT transistor can be function with two different 
regimes (linear regime and saturation regime). 


2.2.1. The linear region (Vg, < Vasat) 

View to variation of potential along the 2-DEG channel and with describe it at each position x of the 
channel, the voltage V,, replace by Vj, — Vc(x) in the expression (1) of sheet carrier concentration of two- 
dimensional electron gas formed at the AlGaN/GaN heterojunction. 


2 


-Kz (r+ [KHAK T m) (Vgs -Ven Tm) -V (x) Ka (T) 


ng(T,m, x) = AG (8) 
Using the following approximation [13]. 
4. K,(T,m).Vo(x) << K3(T) + 4. K4(T,m) (Vys — Ven (T,m) — K, (T)) (9) 


Referring to (8) a simple expression of the concentration n,(T,m). This later independent of the 
position x: 


2 


-kz (+ [KHAK T m) (Vgs -Ven Tm) -K1 (7) 


ns(T,m) = 2.K4(T,m) 


(10) 


Consequently, the drain current from (8) will be expressed by (11). 
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la,(T,m, x) = Z. qu). (nT, m. =Z) (11) 


d. 


We substitute (7) and (10) in the drain current (11), a new equation will appear which describes the drain 
current at each position x of the channel. 


2 
-Ka (T)+ |Z (1) +4.Ka (om). (Vgs -Ven (m) -K1 (T) eta 
las (Tm, x) = Zool 2.Kq(Tam) (aerate) 2) 
Ecvsat ` dx 


The Figure 2 illustrates the equivalent electrical circuit diagram of transistor HEMT AlGaN-GaN, 
where Reis the resistance of the 2-DEG channel, Rg and R, are the parasitic resistances of the drain and 
source, respectively. 


Vys Vas 
EEE aa 


me) 
Ca % 


| 
| 


Oo 
ty 
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Figure 2. Equivalent electrical schema of the HEMT 


From shema in Figure 2 we can conclude the channel potential, for x = 0 and x = L, and it can be given by 
(13) and (14). 
Ve(x)lx=0 = las- Rs (13) 
Ve@)lx=n = Vas — (as + Ra) (14) 


The different step of the integration of the (12), along the channel’s transistor, and with above boundary 
conditions can be given by (15) to (20). 


L (Ho-Ec—Vsat) aVc(x) _ L (adVc(x) 
lush (1 pe Me) dx =2.q.ng(T,m). uo. fy (PE) . ax (15) 
Tas (Dlg + EER.. [Vo(x)]§) = Z. q-ns CT, m). Ho. [Ve C) (16) 


=H, (eiw r, Ẹ ro) Flie (x + HoE Vsa) y+ Z.9.ng(T,m). flo. (Rs + ro) = 


Ec.Vsq Ec-Vsat 
Z.q.Ns(T,m). Ho. Vas = 0 (17) 
We put: 
aq =- (setae cp, FA Ro) (18) 
CYsat 
J ( )) 
—K,(T)+ | K3(T)+4.K4(T).(Vgs—Vin(T,m)-Ky(T) Eo-Vea 
@ = L + Z.q. Ho. (Rs + rol AC ) + tee: t) Vas (19) 
2 
-K2(T)+ le (T)+4.K4(T,m).(Vgs-Vtn(T,m)-K1(T)) 
a3 = —Z.q. to. | ———— (20) 


2.K4(T,m) vids 
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The drain current Ig, expression in linear regime, will be given by the root of (17), and can be illustrated by 
following form (21). 


wat |a3—4.04.03 (21) 


i 
ds 2a, 


2.2.2. The saturation region (Vas > Vasat) 

The free charge carriers (holes and electrons) acquire a velocity of displacement v proportional to 
the applied field, by application an electric field in semiconductor. This velocity can be described by the 
following relation (22). 


We) 


= = Ho 
Sa KONEG) z (14o Ec2sat) e) dx (22) 
EcYsat ` dx 


At high electric field values, the velocity of the carriers is saturated at the value u,,;. Which have allowed us 
to replace the velocity v by Vsat in expression (12) of the drain current Ias. 


2 


-K2(T)+ [E Draka (73) Vgs -Vin Tm- (1) 


2.K4(T,m) 


lasat = Z.q.Vsat (23) 


In the same way, we replace the voltage Vgs by the saturation voltage V4sat in the equations of az 
and az . In this case, the drain current expression in the regime of saturation can be written in the below form. 


—B2+ | B3—-4.04.B3 (24) 


lasat = 2a 
With, 

By = 6, + POEET V4.4 (25) 
and, 

B3 = 63° Vasat (26) 
where, 


2 
-K3 (ryt [KZ +4-Ka(Tm).(Vgs-Ven (T,m)-K1(T)) 


2.K4(T,m) 


62 = L + Z.q. Ho. (Rs + Ra) (27) 


2 
-K3 (T)+ fe (T)+4.K4(T,m).(Vgs-Ven(T,m)-K1(T)) 


2.K4(T,m) 


63 = -Z, q. Ho. 
(28) 


To determine the expression of the saturation voltage V4sat we equalize the (23) and (24) of the 
drain current lgs- 


2 


—Bot [634.01 Bs —K2(T)+ [KZ (1)+4.K4(T,m).(Vgs—Ven(7.m)-Ky (T) 
RENATE ge Z 2 ( ) (29) 


= 24.9.U > 
2.04 q- Vsat 2.K4(T,m) 


We put: 
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2 
-K2(T)+ [Eraka (90) (Vgs—Ven(-)-Ki (7) 


6, =2- d1: Zq Vsat 2KA(Tm) (30) 
In (29) becomes (31). 

—B, + V p2 — 4.1.83 = 6; (31) 
According to the expressions of f, and $3, (31) becomes (32). 

(uoEc—Ysat) (uoEc—vsat) $ 

by = — (5 ee Vac) + a a) — 4.0. a (32) 
As a result, the expression of the saturation voltage Vasat can be expressed by (33). 

Ve -54.(2.5+63) (33) 

sat T 


Ho-Ec-vsat) ) 
2.6 WoFo-vsat) 4 q 6: 
( 1 Ec-Ysat He 


2.3. Low-field mobility uo 

The low field mobility 4ọ has a significant role in the different performances of transistors [16]. 
Based on the simulation of this thermal model, we have observed that the low field mobility varies with 
temperature T, which is presented in the (34). 


Mo(T) = Uo(300) — (1074 x (T — 300)) (34) 


With 449(300) presents the mobility at a temperature equals 300 Kelvin. We deduce that, the low field 


mobility is influenced directly by the temperature, when the temperature! rises, the mobility declines, which 
result to a reduction in the performance of the transistor. 


2.4. The saturation drift velocity v.,, 

Due to the high electric field in electronic devices, the saturation drift velocity Vsat plays an 
important role [17]. According to the (35) below, we have proposed during the simulation of this thermal 
model the variation of the saturation drift velocity v,,;(T) as function of temperature, we can confirm that 
while the temperature Tincreases, the saturation drift velocity decreases, which limits the correct functioning 
of the transistor HEMT in the applications at high temperatures. 


Vsat(T) = Vsar(300) — (10 x (T — 300)) (35) 
With v.q¢(300) presents the saturation drift velocity at to the temperature of 300 Kelvin. 


2.5. Thermal conductivity 

The quantity of heat transferred in a material is known as the thermal conductivity, which can be 
expressed in (W.K"!.m'). The thermal conductivity has modeled by a power law as function of temperature 
as (36). 


300\% 
K(T) = ko (©) [18] (36) 
Where a is the power law coefficient, and kg present the thermal conductivity at 300 Kelvin. In this 
manuscript we will focus on the thermal conductivity of the two substrates: silicon and sapphire, then we 
study the variation of the temperature in the device using these two substrates (silicon and sapphire), and we 
observe the substrate that has a lower increase of the temperature in the transistors HEMTs. For this, we use 


the two following as (37) and (38). 


300 
T 


Keapnire(T) = 0.49 (=Œ) [19] (37) 
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14 
Ks(T) = 1.57 (=) © [20] (38) 
2.6. Self-heating effect 

The power dissipation appears due to the circulation of the drain current in the region of the 2-DEG 
channel of transistor HEMT. This power produces a heating in the 2-DEG channel that is called the self- 
heating that is considered the direct reason for the degradation of the mobility of the carriers and the 
saturation velocity [21]. Consequently, the self-heating is called the channel temperature, its expression can 


be given by the following form [21]. 


1- (2 Pigg" 
i AP / +T (39) 
(Eais m sub sub 
dies 


Tchannel (T,m) = 


Where, P4iss and Toup present successively the power dissipation and the substrate temperature in the 
transistor HEMT, they can be expressed as [21]. 


Paiss(T,m) = Igs(T,m) + Vas (40) 
Tsup(T,m) = 300+ A- Paiss(T,m) (41) 


With, À is a parameter depending on the device thermal resistance [13], and Pp is called the quantity 
with the power dimension (watts) that can be given by (42). 


Po = ere (42) 


in(= su) 


Where, K(T5y), tsup are the thermal conductivity, and the thickness of the substrate respectively, The 
important steps that we will follow to determine the current voltage characteristics considering the effects of 
self-heating are: firstly we calculate the current voltage characteristic (Ias, Vas) in both region (linear and 
saturate) in ambient temperature using (6) and (10). Then we substitute the product of Ig, and Vg, in (40) to 
calculate the power dissipated. Finally, we deduce the channel temperature with which we calculate the 
current-voltage electrical characteristics influenced by the temperature. 


3. RESULTS AND DISCUSSION 


: ; Era 1 
Figure 3 presents, the inverse thermal conductivity Kaa 
sub 


different substrates: sapphire and silicon. According to this figure, we can observe that the inverse of thermal 
conductivity of sapphire increases more than silicon when the temperature T increases. It noticed that thermal 
conductivity of the silicon has less influence with temperature variation than the sapphire substrat. 

The channel temperature variation as a function of the drain voltage Vg, for various values of the 
gate voltage V,, is showing in the Figure 4, from (40), as long as the drain voltage rises, the power dissipation 
increases which translates by an increase of the channel temperature. 

Figure 5 illustrates the current-voltage (Ias, Vas) characteristics of the HEMT transistor for various 
variables of gate voltage V,, from 1.5 to -2.5V with a step of 2V using the sapphire substrate and taking in 
consideration the existence of self heating and the absence of this latter. From this figure, we note that the 
experimental results [22] shown through the violet dots agree well with our theoretical calculation (in 
presence of self-heating) confirming the validity of our model. From the same figure we observe a 
discrepancy AJ between the current-voltage electrical characteristics with and without self-heating. This can 
be explained by the growth of the 2-DEG channel temperature, this latter confirmed well by (40) when the 
drain voltage Vg, increases, which allows the discrepancy AJ appearances. 

Figure 6 illustrats the self-heating effects developed by the variation of the drain current Ig, as a 
function of the drain voltage Vas while the value of the gate voltage V,,equaled zero V,, = 0 for the substrats 
sapphire and silicon in the HEMT transistor AlGaN-GaN. It is observed that the production of the devices 
with the silicon substrate has a high current density instead of using sapphire substrate. This indicate the 
importance of the thermal management to reach the higher performances of the devices. The calculated 
results shows a good agreement with the experimental results [23], which means the validation of the 
proposed model. 


as a function of the temperature T for 
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In the present calculation the different parametrs used to simulate this current developed model are 
illustrated in the Table 2. Also the values of parameters K,, K, and K3 at different temperatures T which 
calculated by the same method [24]. Using the effective mass m* = 0.22mg [25], used too in this developed 


model are obtained in Table 3. 


—e— Sapphire 
—<— Silicon 


1/(K(T.,,) (CMK/W)) 


0,0 L L a L a l 4 
100 200 300 400 500 600 


Temperature T(K) 


Figure 3. The evolution of the reverse thermal 


(cmK/W)) as a function of the 
K(Tsub) 


temperature T for different substrates silicon and 
sapphire 


conductivity 
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Figure 5. The drain current Ig, variation versus the 
drain voltage V,, for various values of the gate 
voltage V,,, taking in consideration the absence and 
the presence of self-heating 
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Figure 4. The channel temperature variation versus 
the drain voltage Vgs for various values of the gate 
voltage V,, for transistor HEMT 
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Figure 6. The drain current Jg, variation versus the 
drain voltage V,, for sapphire and silicon substrate 
materials for L = 0.45 um and m = 0.2 


Table 2. Different parameters used in the present simulation 


Parameter Description Figure5 Figure 6 
dep(nm) Thickness of the undoped cover layer 3 4 
d,(nm) Thickness of the doped layer 15 16 
dsp(nm) Thickness of the undoped layer 7 6 
Z(um) Gate width 15 50 
L(um) Gate length 1.5 0.45 

R,(Q) Parasitic source resistance 0.6 0.6 
Ra(Q) Parasitic drain resistance 0.9 0.9 
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Table 3. Values of the parameters K,, Kz and K3, which calculated by the same method [24] using the 
effective mass m* = 0.22m [25] 


T(K) KV) K,(Vm) K,(Vm*) 
300 -0.2775 1.68742x 10 -4.2371x108 
350 -0.328 1.9083x10* -5.2147x10"8 
400 -0.3792 2.1232x10% -6.1660x10% 
450 -0.4313 2.3357x10% -7.1058x10"8 
470 -0.4512 2.41030x10ë -7.4352x10!8 


4. CONCLUSION 

In this work, we investigated the effect of the self-heating on current-voltage electrical 
characteristics for the transistor HEMT based on GaN. Firstly, we studied the self-heating effect for one 
substrate, then, we used two substrates to choose the coldest substrate to be used in a high voltage 
functionning. To highlight the validity and the confirmation of our proposed model, we conducted a 
comparative study between the simulation results and experimental measurements, this study proved that 
there is a good match between them, which confirms the validity of our proposed model. 
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